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Final  report:  Scalable  directed  self-assembly  using  ultrasound  waves 

Army  Research  Office  STIR  contract  W91  INF 14 10565 
Bart  Raeymaekers,  Assistant  Professor 
Department  of  Mechanical  Engineering  -  University  of  Utah 

Synopsis 

We  aim  to  understand  how  ultrasound  waves  can  be  used  to  create  organized  patterns  of 
nanoparticles  in  a  host  medium  such  as  a  polymer  matrix  material.  The  critical  difference  between 
the  ultrasound  technology  studied  in  this  project,  and  other  directed  self-assembly  techniques  is  its 
scalability  and  flexibility.  During  this  9  month  project,  we  have  derived  the  theory  and 
implemented  a  model  to  simulate  how  patterns  of  nanoparticles  dispersed  in  a  host  medium  are 
assembled  by  means  of  standing  ultrasound  waves.  Additionally,  we  have  obtained  experimental 
validation  of  this  model.  This  research  finds  application  in  design  and  processing  of  novel  bulk 
materials  and  could  lead  to  processing  of  bulk  nanocomposite  materials  with  tailored  properties. 

Contribution  to  the  Army’s  mission  and  benefit  to  society 

The  ability  to  create  user- specified,  organized  patterns  of  nanoparticles  in  a  macroscale  volume 
of  polymer  material  is  anticipated  to  expedite  the  future  bulk  processing  of  polymer-based  multi¬ 
functional  nanocomposite  materials  with  tailored  properties.  Nanocomposite  materials  attempt  to 
take  advantage  of  the  unique  physical  and  mechanical  properties  of  nanoparticles  (e.g.  carbon 
nanotubes)  to  create  ultra-high  strength-to-weight  and  stiffness-to-weight  ratio  materials  with 
mechanical  properties  that  outperform  those  of  state-of-the-art  materials,  at  a  fraction  of  the  cost. 
These  materials  could  be  used  to  replace  traditional  materials  such  as  metals,  to  reduce  strain  on 
the  Warfighter  by  reducing  payload.  Lightweight  materials  could  furthermore  improve  energy 
efficiency  of  rolling  stock,  to  sustain  dominance  in  military  operations  over  land.  In  addition, 
society  could  benefit  from  these  strong,  lightweight  materials  in  the  context  of  energy  efficiency 
of  vehicles,  amongst  other  civilian  applications.  In  addition,  the  ability  to  tailor  patterns  of 
nanoparticles  embedded  in  a  polymer  matrix  by  means  of  ultrasound  directed  self-assembly,  opens 
the  door  to  designing  and  manufacturing  multi-functional  materials  with  tailored  properties, 
beyond  improving  mechanical  properties  and  manufacturing  lightweight  materials. 

1.  Statement  of  the  problem  studied 

The  scientific  objective  of  this  research  is  to  understand  how  ultrasound  waves  can  be  used  to 
create  large-scale,  user-specified  patterns  of  nanoparticles  dispersed  in  a  host  medium.  The 
overarching  engineering  objective,  to  which  this  research  applies,  is  to  manufacture  macroscale 
polymer-based,  multi-functional  nanocomposite  materials  with  tailored  properties. 

We  have  formulated  the  following  two  research  tasks. 

Task  1:  Derive  the  theory  and  implement  a  model  to  simulate  organizing  nanoparticles  into  user- 
defined  orientations  and  patterns  in  a  host  medium  using  ultrasound  directed  self-assembly  (DSA). 
Anticipated  result  and  deliverable:  A  scientific  tool  to  predict  particle  organization,  pattern,  and 
orientation,  based  on  the  operating  and  design  parameters  of  the  ultrasound  transducers,  and  the 
properties  of  the  particles  and  host  medium  material. 

Task  2:  Experimentally  demonstrate  that  user-specified  patterns  and  orientations  of  nanoparticles 
dispersed  in  a  host  medium  can  be  obtained  using  ultrasound  directed  self-assembly. 

Anticipated  result  and  deliverable:  Validation  of  the  theory  and  model  developed  in  Task  1. 
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2.  Summary  of  most  important  results 

During  the  9  month  STIR  contract,  we  have  achieved  the  following  key  results: 

•  We  have  derived  the  theory  and  implemented  a  model  to  simulate  the  physics  of  two- 
dimensional  ultrasound  directed  self-assembly. 

•  We  have  packaged  the  ultrasound  directed  self-assembly  model  into  a  scientific  tool  with 
user-friendly  graphical  user  interface  (GUI).  The  GUI  allows  the  user  to  solve  any 
ultrasound  directed  self-assembly  problem  in  two  dimensions,  without  having  to 
understand  the  underlying  model. 

•  We  have  experimentally  validated  our  ultrasound  directed  self-assembly  model  in  water 
and  polymer. 

In  the  remainder  of  Section  2,  we  will  discuss  each  of  the  key  results  in  greater  detail. 

2.1  Developing  an  ultrasound  directed  self-assembly  model 

Enabling  the  practical  implementation  of  ultrasound  DSA  as  a  manufacturing  technique 
requires  linking  the  transducer  arrangement  and  settings  that  generate  the  standing  ultrasound 
wave  field,  to  the  resulting  patterns  of  nanoparticles.  This  translates  into  two  problems:  (1)  the 
forward  problem  entails  calculating  the  pattern  of  nanoparticles  resulting  from  user-specified 
transducer  settings,  and  (2)  the  inverse  problem  involves  calculating  the  transducer  settings 
required  to  achieve  a  user-specified  pattern.  In  both  cases  the  transducer  arrangement  is  also  user- 
specified. 

The  forward  problem  is  solved  by  computing  the  acoustic  radiation  force  associated  with  the 
standing  ultrasound  wave  field.  The  resulting  pattern  of  nanoparticles  is  then  calculated  as  the 
stable  fixed  points  x/  of  the  acoustic  radiation  force,  defined  as  the  location(s)  where  the  force  is 
zero  and  points  toward  x/  in  the  surrounding  region  [1].  The  inverse  problem  is  solved  either 
directly  or  indirectly,  for  specific  reservoir  geometries,  transducer  arrangements,  and  patterns  of 
nanoparticles.  For  circular  reservoirs  lined  with  transducers  that  are  equipped  with  backing  and 
absorbing  layers  to  remove  reflected  ultrasound  waves,  the  transducer  phases  to  create  dot  patterns 
are  calculated  directly  as  a  function  of  the  desired  geometry  [2,3].  Alternatively,  the  inverse 
problem  is  solved  indirectly  by  first  solving  the  forward  problem  for  a  range  of  possible  transducer 
settings  to  create  a  “map”  of  achievable  patterns  as  a  function  of  the  transducer  settings  [4-6]. 
However,  these  existing  methods  fail  to  provide  a  general  solution  to  the  inverse  problem,  because 
a  unique  function  or  map  must  be  derived  for  each  specific  reservoir  geometry,  transducer 
arrangement,  or  user-specified  pattern  geometry.  Conversely,  using  ultrasound  DSA  as  a 
manufacturing  technique  requires  the  ability  to  calculate  transducer  settings  for  any  feasible  user- 
specified  pattern  of  nanoparticles,  using  a  reservoir  with  an  arbitrary  geometry  and  transducer 
arrangement.  Theoretically,  an  infinite  number  of  transducers  enables  assembly  of  any  user- 
specified  pattern.  However,  fewer  transducers  limit  the  patterns  that  are  feasible  to  assemble  with 
the  ultrasound  wave  field. 

Thus,  we  demonstrate  a  direct  solution  methodology  to  the  inverse  problem  for  any  feasible 
user-specified  pattern  within  a  2D  reservoir  with  an  arbitrary  simply-closed  geometry  and 
transducer  arrangement.  We  accomplish  this  by  linking  the  transducer  settings  to  the  obtained 
patterns  in  two  steps.  First  we  calculate  the  standing  ultrasound  wave  field  in  an  arbitrary  shaped 
reservoir  lined  with  transducers  around  its  perimeter  as  a  function  of  the  transducer  settings,  using 
the  boundary  element  method  based  on  Green’s  third  identity  [7],  which  relates  the  wave  field 
within  a  simply-closed  domain  to  the  boundary  conditions  imposed  on  the  perimeter  of  the  domain. 
Then,  we  apply  Gor’kov’s  method  [1,8]  of  calculating  the  acoustic  radiation  force  acting  on  a 
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spherical  particle  to  determine  the  pattern  of  nanoparticles  resulting  from  the  standing  ultrasound 
wave  field.  The  inverse  problem  is  then  posed  as  an  optimization  problem  whose  solution  is  found 
using  eigendecomposition.  A  theoretical  derivation  and  experimental  validation  is  presented. 

Figure  1  shows  a  2D  reservoir  of  arbitrary  geometry  filled  with  a  fluid  medium  of  density  pm 
and  sound  speed  Cm,  and  with  Nt  transducers  of  density  pt  and  acoustic  impedance  Zi  around  the 
perimeter.  The  inset  of  Fig.  1  shows  a  magnified  view  of  a  section  of  the  reservoir,  illustrating  the 
discretization  of  the  boundary  into  Nb  boundary  elements  {Nb  >  Nt)  and  the  domain  into  Nd  domain 
points  (red  dots),  which  can  be  selected  in  any  arrangement.  The  j*  boundary  element,  identified 
by  its  center  point  q/  is  e(q/)  wide,  and  is  driven  by  the  transducer  setting  v(q/),  i.e.,  the  complex 
harmonic  velocity  (amplitude  and  phase)  of  the  transducer  surface  along  the  normal  direction  n(q/). 
Additionally,  a  test  point  x/  within  the  reservoir  is  indicated  with  respect  to  the  reservoir  origin  o. 

We  use  the  boundary  element  method  to 
calculate  the  ultrasound  wave  field  with 
frequency  coo  in  terms  of  the  time- 
independent,  complex  scalar  velocity 
potential  cp,  at  each  domain  point  within  D. 

We  note  that:  (1)  ^  must  be  a  wave  and,  thus, 
satisfy  the  Helmholtz  equation  (y^cp  +  k^^cp  = 

0)  in  D,  where  ko  =  coolcm  is  the  wave  number 
of  the  ultrasound  wave  field  in  the  fluid- 
medium.  (2)  The  impedance  boundary 
condition  dcpidn  +  ikoZcp  =  iptlpm)v  must  be 
satisfied  on  S,  where  Z  =  pmcJZt  is  the 
impedance  ratio,  accounting  for  the 
absorption  and  reflection  of  the  ultrasound 
wave  within  the  fluid-medium  as  it  interacts 
with  the  transducers.  (3)  cp,  ^cp,  and  ^^cp  must 
be  continuous  in  D  and  on  5,  except  where 
the  boundary  is  non-smooth  (e.g.  sharp 
corners  of  the  reservoir),  satisfying  the 
assumptions  of  Green’s  third  identity  [7].  By 
arranging  the  values  of  v  for  each  transducer  into  a  vector  v  the  ultrasound  wave  field  in  D  is 
calculated  as  [9] 

(p  =  PWv.  (1) 

The  matrix  W  maps  each  boundary  element  to  its  corresponding  transducer,  i.e.,  Wjm  =  1  if  the  f' 
boundary  element  is  contained  within  the  m*  transducer,  otherwise  wjm  =  0.  Additionally,  each 
term  pij  of  the  matrix  P  corresponds  to  the  ultrasound  wave  field  created  in  x/  by  a  point  source 
located  at  q,  on  5,  including  all  reflections  from  the  reservoir  walls.  All  pij  terms  are  calculated  in 
matrix  form  as 


Figure  1 :  2D  fluid-filled  reservoir  with  arbitrary  geometry,  lined 
with  Nt  transducers.  The  insert  illustrates  the  discretization 
scheme  used  by  the  boundary  element  method,  dividing  the 
boundary  into  Nb  boundary  elements  (black  dots)  and  the  domain 
into  Nd  domain  points  (red  dots).  Additionally,  the  insert  shows 
the  length  e(q/)  and  normal  direction  n(q;)  of  the  boundary 
element  centered  at  qj. 


F  = 


B-A(4I  +  A)  ‘b 


(2) 


for  all  Nd  domain  points.  Here,  I  is  the  identity  matrix  and  each  term  aij  and  bij  of  the  matrices  A 
and  B  are  calculated  as 
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Pm 

i  =  (-1)^^^  is  the  unit  imaginary  number,  the  function  5(q,,  x/)  =  0  when  q,  =  x/,  otherwise  it  is  1,  and 
G(q/,x/)  is  the  Green’s  function,  which  represents  the  free-field  ultrasound  wave  emitted  from  a 
point-source  located  at  q/  and  measured  at  location  xi ,  defined  in  2D  as 

G(q;,X;)  =  -^//o(A:o|qy-x,|).  (5) 

Here,  Ho  is  the  0*  order  Hankel  function  of  the  first  kind,  and  Iq,  -  x/l  is  the  Euclidean  distance 
between  points  q,  and  xi.  The  matrices  A  and  B  in  Eq.  (2)  are  calculated  using  the  same  formulas 
as  A  and  B,  differing  only  by  the  selection  of  the  points  x/,  which  lay  on  S  for  A  and  B,  and  lay  in 
D  for  A  and  B.  Thus,  using  the  boundary  element  method  we  link  the  transducer  settings  to  the 
resulting  ultrasound  wave  field. 

To  link  the  ultrasound  wave  field  to  the  final  pattern  of  nanoparticles,  we  calculate  the  acoustic 
radiation  force  acting  on  a  nanoparticle  of  density  pp  and  sound  speed  Cp,  dispersed  in  a  fluid 
medium  at  location  xi  in  D,  as 

f/=-vc/,,  (6) 

where  Ui  is  the  acoustic  radiation  potential  in  x/.  Eor  a  spherical  particle  with  radius  Vp  «  Xo,  where 
/lo=  Inoiolcf,  Ui  is  calculated  as  [1,8] 

C/,=v"Q,y,  (7) 

where  v^  is  the  conjugate  transpose  of  v,  and  the  matrix  Q/  is  calculated  as 


Here,  p/^  is  the  row  of  P,  and  =  HpmCm  and  Pp  =  MppCp  is  the  compressibility  of  the  fluid 
medium  and  particle,  respectively.  From  Eq.  (7),  the  final  pattern  of  nanoparticles  is  calculated  as 
the  region(s)  consisting  of  points  x/,  where  Ui  is  locally  minimum.  Thus,  to  achieve  a  user-defined 
pattern  consisting  of  the  set  of  desired  points  Xides,  each  value  Ui,  corresponding  to  each  point  x/  G 
'Xides,  must  be  locally  minimum.  To  pose  this  as  an  optimization  problem  with  a  single  objective 
function,  we  relax  the  requirement  of  local  minimality.  Instead,  we  aim  to  minimize  the  average 
value  of  Ui  over  all  points  x/  G  Xides,  which  is  written  as  the  quadratic  functions 

(j  =  y^Qy,  (9) 

/V 

where  the  matrix  Q  is  the  average  of  the  matrices  Q/  corresponding  to  each  desired  point  x/  G  Xdes- 
The  matrix  Q  is  indefinite  because  [7  has  no  lower  bound.  Physically,  this  occurs  because  particles 
can  be  trapped  at  the  desired  positions  more  effectively  by  increasing  the  amplitudes  of  the 
transducers  indefinitely  (lv|  ^  co).  Practically,  the  function  generator  used  to  excite  the  transducers 
limits  the  transducer  amplitudes  to  finite  values  and,  thus,  we  constrain  the  magnitude  Ivl  =  a, 
where  a  is  a  real,  scalar  value  representing  the  maximum  harmonic  velocity  that  can  be  achieved 
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with  a  function  generator  and  transducer.  Hence,  we  formulate  the  constrained  quadratic 
optimization  problem 

min  U  ,  subject  to  |v|  =  a .  (10) 

From  Eq.  (10),  the  optimal  transducer  settings  v*  are  calculated  as  the  eigenvector  of  length  a 
corresponding  to  the  smallest  eigenvalue  of  Q.  Thus,  we  solve  the  inverse  problem,  finding  the 
transducer  settings  v*  required  to  achieve  a  feasible  user-specified  pattern  of  nanoparticles,  defined 
as  a  set  of  desired  points  Xdes. 

2.2  Packaging  the  ultrasound  directed  self-assembly  model  as  a  scientific  tool 

We  have  packaged  this  model  in  a  user-friendly  graphical  user  interface  (GUI).  The  GUI 
enables  solving  both  the  forward  and  inverse  problem,  and  simulating  ultrasound  directed  self- 
assembly,  without  the  need  to  handle  the  source  code  of  the  model.  The  GUI  allows  selecting  the 
reservoir  dimensions,  transducer  arrangement  and  properties,  host  medium  and  particle  properties. 
In  the  case  of  solving  the  forward  problem,  the  model  calculates  the  pattern  of  nanoparticles  that 
is  obtained  for  the  selected  properties,  and  for  user-specified  transducer  inputs.  Conversely,  in  the 
case  of  solving  the  inverse  problem,  the  model  finds  the  transducer  settings  needed  to  obtain  a 
user-specified  pattern  of  nanoparticles  for  the  selected  properties.  Figure  2  shows  a  screenshot  of 
the  scientific  tool  and  GUI  that  we  have  developed  to  package  the  ultrasound  DSA  model  and 
enable  easy  access  to  running  simulations.  We  have  included  two  movies  illustrating  the  use  of 
the  GUI  to  solve  the  forward  problem  (Appendix  1)  and  the  inverse  problem  (Appendix  2). 
In  addition,  we  have  included  a  ZIP  file  with  all  source  code  and  an  executable  file  that  allows 
running  the  scientific  tool  across  different  platforms  (Appendix  3). 


ll  MainWindow 


I 


Visualization 
window 

•  Sketch 
desired 
pattern 

•  Observe 
results 


[  Save  image  | 


I a  T 


Mode  selection 

(forward/inverse 

problem) 

Post-processing 
Input  parameters 


Figure  2:  Scientific  tool  developed  to  package  the  ultrasound  directed  self-assembly  model,  showing  the  graphical  user  interface 
(GUI)  that  enables  solving  the  forward  and  inverse  problem,  and  post-process  the  results. 
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Figure  3:  Experimental  apparatus.  A  user-specified  pattern  is 
defined  within  the  modeled  reservoir,  then  the  inverse  problem 
is  solved,  yielding  the  optimal  transducer  settings  which  are 
applied  to  the  machined  reservoir,  generating  the  actual  pattern 
of  nanoparticles.  A  CCD  camera  is  then  used  to  image  the  actual 
pattern. 


2.3  Experimental  validation  in  water 

Figure  3  shows  a  schematic  of  the 
experimental  apparatus  that  we  have  used  to 
validate  our  model.  A  12.75x12.75  mm 
square  reservoir  containing  water  with 
dispersed  80  nm  carbon  nanoparticles  is 
simulated  using  our  model  described  in 
Section  2.1,  and  manufactured  to 
experimentally  validate  the  model.  In  both 
the  model  and  experiment,  each  of  the  four 
reservoir  walls  consists  of  a  single  PZT 
transducer  (type  SMI  12)  with  center 
frequency  1.5  MHz,  driven  by  the  complex 
transducer  settings  (amplitude  and  phase)  vi, 

V2,  V3,  and  V4,  respectively.  The  optimal 
transducer  settings  needed  to  create  a  feasible 
user-specified  pattern  are  computed  by 
solving  Eq.  (10)  and  then  applied  to  the 
transducers  in  the  experimental  setup  to 
validate  the  actual  pattern  of  nanoparticles. 

Figure  4  shows  the  two  types  of  patterns 
that  can  be  achieved  in  a  square  reservoir  with  one  transducer  per  side.  These  patterns  consist 
either  of  (a)  lines  spaced  XqH  apart,  oriented  along  either  the  x-  or  y-directions,  or  (b)  dots  arranged 
in  a  square  grid  formation  spaced  XqH  apart  in  both  the  x-  and  y-directions.  The  user-defined  pattern 
is  depicted  in  red,  from  which  the  complex  transducer  setting  (amplitude  and  phase)  are  calculated 
using  the  “inverse”  method  explained  in  section  2. 1  of  this  report.  These  settings  are  then  applied 
to  the  respective  transducers,  and  the  80 
nm  carbon  nanoparticles  assemble  into 
the  user-specified  pattern.  A  picture  of 
the  experimentally  obtained  pattern  is 
shown,  and  an  inset  shows  a  magnified 
view  of  the  experimental  results, 
overlaid  with  the  user-specified  pattern. 

The  line  and  dot  patterns  can  also  be 
shifted  in  the  x-  or  y-direction.  To 
validate  our  method  of  solving  the 
inverse  problem  in  an  unconstrained 
fashion,  we  have  prescribed  line  and  dot 
patterns  at  positions  shifted  along  the  x- 
and  y-  direction  over  a  distance  of  Ax  G 
[0,  7o/2),  demonstrating  that  user- 
specified  patterns  can  be  created 
throughout  the  reservoir.  For  each 

pattern  we  compare  the  user-specified  Figure  4:  Experimental  results  indicating  the  desired  patterns 
pattern  to  the  experimentally  obtained 

^  r  j  transducer  settings  calculated  using  Eq.  (9)  lor  a  desired  (a)  line 

pattern  by  calculating  the  pattern  error  pattern  and  (b)  dot  pattern  of  nanoparticles. 
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Figure  5:  Pattern  error  for  desired  line  and  dot  patterns  as  they  are  shifted  the  x-  and  y-directions  in  increments  of  Ax/An  =  0.0625. 
Additionally,  images  of  the  desired  dot  and  line  patterns  are  superimposed  over  the  actual  dot  and  line  patterns  for  the  shifted 
distances  Ax/Xo  =  {0.125,  0.250,  0.375}. 


Epat  as  the  average  distance  between  the  centers  of  the  user-specified  and  experimentally  obtained 
pattern  features  (lines  or  dots),  normalized  by  the  nominal  pattern  spacing  Xo/2.  Figure  5  shows  the 
pattern  error  as  a  function  of  the  normalized  pattern  shift  distance  Ax/Xo  for  line  (triangle  marker) 
and  dot  (dot  marker)  patterns.  Additionally,  Fig.  5  shows  images  of  the  user-specified  patterns 
superimposed  on  the  experimentally  obtained  patterns,  for  line  and  dot  patterns  that  have  been 
shifted  by  =  {0.125,  0.250,  0.375}  in  the  x-  and  y-  directions.  We  observe  that  the  pattern 
error  Epat  remains  less  than  16.0%  and  17%  for  line  and  dot  patterns,  respectively,  indicating  good 
agreement  between  the  user-specified  and  experimentally  obtained  patterns.  In  practice,  these 
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patterns  are  not  identical  due  to  imperfections  in  the  experiment  and  assumptions  in  the  model. 
The  reservoir  is  not  perfectly  square,  and  the  transducers  may  be  slightly  misaligned,  which  affects 
the  standing  ultrasound  wave  field.  Additionally,  the  model  neglects  viscous  and  thermal  effects, 
which  can  produce  acoustic  streaming  and  alter  the  fluid  properties,  which  are  not  accounted  for 
in  the  model. 

When  adding  additional  transducers  to  the  setup,  more  exotic  patterns  can  be  created.  Figure 
6  depicts  two  examples  of  patterns  that  are  obtained  with  two  independent  transducers  per  side  of 
a  square  reservoir,  i.e.,  8  transducers  in  total.  The  user-defined  pattern  is  depicted  in  red,  from 
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Figure  6:  Exotic  patterns  obtained  with  8  transducers,  (a)  combined  line  and  dot  patterns  and  (b)  a  shifted  line  pattern.  Red  line 
depicts  user-specified  pattern,  and  picture  shows  experimental  implementation. 

which  the  complex  transducer  setting  (amplitude  and  phase)  are  calculated  using  the  “inverse” 
method  explained  in  section  2.1  of  this  report.  These  settings  are  then  applied  to  the  respective 
transducers,  and  the  80  nm  carbon  nanoparticles  assemble  into  the  user-specified  pattern.  A  picture 
of  the  experimentally  obtained  pattern  is  shown,  and  an  inset  shows  a  magnified  view  of  the 
experimental  results,  overlaid  with  the  user-specified  pattern. 
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2.4  Experimental  validation  in  polymer 

We  have  used  ultrasound  directed  self-assembly  to  manufacture  dogbone  specimens  of 
polymer  matrix  composite  materials  with  aligned  carbon  nanotubes  (CNTs)  as  reinforcement.  We 
have  focused  on  manufacturing  materials  with  a  CNT  weight  fraction  in  excess  of  1%,  as  it  is  well- 
known  that  small  weight  fractions  have  a  limited  effect  on  the  mechanical  properties  of  the 
material.  A  dogbone  shaped  reservoir  is  machined  from  high  density  polyethylene,  which  displays 
low  adhesion  to  the  polymer  matrix  material  and  allows  easy  extraction  of  the  finished  composite 
material  specimens  from  the  reservoir.  The  reservoir  is  mounted  on  a  removable  glass  base.  Two 
parallel  piezoelectric  transducers  made  of  PZT-4,  with  resonant  frequency  of  1.5  MHz,  are 
embedded  in  recess  slots  cut  in  the  walls  of  the  reservoir,  and  oppose  each  other  over  the  entire  18 
mm  gauge  length.  The  resulting  dogbone  specimen  resembles  an  ASTM  D638  Type  V  plastic 
sample,  and  is  sufficiently  large  (70  mm  length)  to  illustrate  that  macroscale  components  can 
successfully  be  manufactured  using  the  ultrasound  alignment  method.  Multi-walled  carbon 
nanotubes  (MWCNTs)  as  opposed  to  single-walled  carbon  nanotubes  (SWCNTs)  are  used  in  the 
experiments  because  of  cost.  Figure  7  illustrates  the  manufacturing  and  alignment  process. 
MWCNTs  with  a  diameter  of  50-80  nm  and  length  of  10-20  pm  are  added  to  part  A  of  a  two-part, 
low-viscosity,  fast  curing  thermoset  resin  (Smooth-Cast  300  thermoset,  Smooth-On  Inc.)  (Fig.  7 
a)),  and  dispersed  by  mwcnts  Pan  a  resin 

means  of  bath  sonication 
(Fig.  7  b)).  Part  B  of  the 
resin  is  then  added  (Fig.  7 

c) )  prior  to  transferring  the 
mixture  to  the  dogbone 
shaped  reservoir  (Fig.  7 

d) ).  The  viscosity  of  the 
resin  at  room  temperature 
is  80  cPs.  A  standing  bulk 
acoustic  wave  (BAW)  is 
created  between  the  two 

opposing  PZT  plates  to 
align  the  MWCNTs  at  the  Figure  7:  Manufacturing  process  of  a  dog-bone  shaped  composite  material  specimen  with 
^  aligned  MWCNTs 

nodes  of  the  standing 

wave  pattern  (Fig.  7  e)),  when  this  force  exceeds  the  drag  force  acting  on  the  MWCNTs.  The  PZT 
plates  are  driven  by  a  function  generator  (Tektronix,  AFG  3102),  amplified  by  a  45  dB  50  W  RF 
power  amplifier  (Electronic  Navigation  Industries,  440LA).  Cross-linking  of  the  mixture  starts 
immediately  following  the  addition  of  the  second  resin  part  to  the  mixture.  The  sample  is  extracted 
from  the  reservoir  upon  completion  of  the  cross-linking  process  (Fig.  7  f)). 

Figure  8  a)  displays  a  typical  example  of  a  dog-bone  specimen  with  one  weight  percent  (wt  %) 
aligned  MWCNTs,  fabricated  using  ultrasound  waves  with  a  frequency  of  1.477  MHz.  Figure  8  b) 
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shows  a  detail  of  the  gauge  section  with  aligned 
MWCNTs,  obtained  with  an  optical  microscope, 
after  locally  polishing  the  surface  of  the  specimen, 
to  reveal  the  presence  of  the  aligned  MWCNTs. 

The  spacing  between  two  rows  of  aligned 
MWCNTs  is  458  pm,  which  is  expected  based  on 
the  sound  speed  of  the  resin,  and  the  occurrence  of 
two  nodes  per  wave  length. 

We  have  performed  quasi-static  tensile  testing 
for  specimens  with  MWCNT  weight  fraction 
between  1  and  10%  and  compared  the  Young’s 
modulus  and  ultimate  tensile  strength  to  virgin  polymer  material,  and  to  polymer  material  with  1 
weight  percent  randomly  oriented  MWCNTs.  We  observed  that  the  ultimate  tensile  strength  and 
Young’s  modulus  increase  with  increasing  CNT  weight  fraction,  and  observed  that  randomly 
oriented  MWCNTs  do  not  significantly  improve  the  mechanical  properties  of  the  polymer  matrix 
material,  in  line  with  previous  research.  However,  we  obtained  a  200%  increase  in  the  ultimate 
tensile  strength  for  the  case  of  10%  weight  fraction  of  MWCNTs  aligned  in  the  polymer  matrix. 


Locally  polished 
area 

Figure  8:  a)  Aligning  MWCNTs  in  a  70  mm  long  dog 
bone  sample,  showing  the  scalability  of  the  acoustic 
alignment  technique  to  align  MWCNTs  in  a  polymer 
matrix  material.  This  sample  contains  5  Vol%  (1  wt%) 
MWCNTs.  A  small  section  is  polished  to  reveal  the 
presence  of  the  aligned  MWCNTs  and  b)  close-up  of 
aligned  MWCNTs 


3.  Conclusion 

We  have  derived  a  method  of  ultrasound  DSA  which,  for  the  first  time,  provides  a  direct 
solution  to  the  inverse  problem  for  a  reservoir  with  arbitrary  geometry,  calculating  the  transducer 
settings  necessary  to  achieve  any  feasible  user-specified  pattern  of  nanoparticles.  In  contrast  to 
existing  methods  in  the  literature,  which  are  unique  to  each  specific  reservoir  and  pattern,  this  new 
method  provides  a  general  technique  for  solving  the  inverse  problem  in  reservoirs  with  for 
arbitrary  geometries,  transducer  arrangements,  and  patterns.  Using  a  square  reservoir,  feasible  line 
and  dot  patterns  are  specified,  and  our  new  method  is  used  to  calculate  the  transducer  settings 
necessary  to  generate  the  patterns.  These  transducer  settings  are  applied  to  an  experimental 
reservoir,  generating  the  actual  patterns  of  nanoparticles,  which  are  then  compared  to  the  desired 
patterns.  Good  agreement  is  found  between  the  desired  patterns  and  the  actual  patterns, 
demonstrating  our  new  method’s  ability  to  effectively  calculate  the  transducer  settings  necessary 
to  generate  user-specified  patterns. 


4.  Research  products 

Scientific  tool:  We  have  included  the  scientific  tool  developed  to  simulate  ultrasound  directed 
self-assembly  in  Appendix  3.  All  source  code  has  been  included,  and  an  executable  file  that 
launches  the  GUI  of  the  scientific  tool. 

Journal  publications:  We  have  focused  during  the  9  month  project  on  completing  the  research 
and  implementing  the  scientific  tool.  While  we  already  have  published  one  journal  publication,  we 
are  still  in  the  process  of  publishing  the  research  results  of  this  project.  We  have  included  a  list  of 
the  different  journal  publications  that  we  have  published  or  are  preparing.  We  include  the  ARO 
funding  in  the  acknowledgments  of  all  these  publications. 


Corbitt  SJ,  Petersen  SJ,  Francoeur  M,  Raeymaekers  B,  State-of-the-art  fabrication  of  Mie 
resonance-based  dielectric  metamaterials  operating  at  optical  frequencies  for  use  in  engineering 
devices:  A  review;  Journal  of  Quantitative  Spectroscopy  and  Radiative  Transfer,  Vol.  229(4),  pp. 
547-556  (2015) 
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Greenhall  JJ,  Guevara  Vazquez  F,  Raeymaekers  B,  Ultrasound  directed  self-assembly  of  2D 
arbitrary  patterns  of  nanoparticles,  in  preparation  for  submission  (2015) 

Greenhall  JJ,  Guevara  Vazquez  F,  Raeymaekers  B,  Two-dimensional  constrained  optimization 
methods  for  ultrasound  directed  self-assembly,  in  preparation  for  submission  (2015) 

Homel  LJ,  Greenhall  JJ,  Raeymaekers  B,  Aligning  ultra-high  weight  fraction  carbon  nanotubes  in 
a  polymer  matrix  using  bulk  ultrasound  waves,  in  preparation  for  submission  (2015) 

5.  Interactions  with  Army  Research  Labs 

LJ  Holmes:  I  have  visited  LJ  Holmes’s  research  group  (Micro-Compositronics  and  Rapid 
Operations  Lab)  in  July  2013  at  Aberdeen  Proving  Grounds  (APG),  to  discuss  a  possible 
collaboration.  The  idea  is  to  integrate  the  ultrasound  directed  self-assembly  technique  with  additive 
manufacturing.  LJ  is  an  expert  in  additive  manufacturing  and  showed  me  around  his  laboratory 
and  several  other  labs  at  APG,  I  also  gave  a  seminar  at  APG  during  that  visit.  During  the  STIR 
grant,  I  have  kept  regular  contact  with  LJ  and  have  solicited  his  feedback  to  ensure  that  our  research 
is  in  line  with  the  Army  mission. 

Steve  Taulbee:  Steve  Taulbee  visited  the  University  of  Utah  in  January  2015  and  I  briefly  met 
with  him  to  discuss  the  Army  Research  on  which  we  are  working.  Steve  invited  me  to  visit  APG 
again,  and  I  have  given  a  seminar  at  APG  in  April  2015,  in  addition  to  meeting  again  with  LJ 
Holmes,  Ray  Brennan,  and  others,  to  discuss  our  research  progress  in  light  of  Army  needs  and 
future  directions. 

Raymond  Brennan:  I  have  had  several  phone  conversations  about  the  fit  and  a  possible 
collaboration  with  the  Energy  Coupled  to  Matter  program,  with  Ray  Brennan.  We  plan  to  extend 
the  research  performed  during  this  STIR  contract  to  use  ultrasound  directed  self-assembly  for 
manufacturing  of  multi-functional  materials,  and  we  will  collaborate  with  Ray  on  this  effort. 
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Appendices 

We  have  provided  Dropbox  links  to  each  of  the  appendices,  as  the  ARO  reporting  website  does 
not  allow  for  submitting  these  extra  files.  However,  we  believe  they  are  an  intricate  part  of  our 
reporting. 

Appendix  1 :  Movie  demonstrating  solving  the  forward  problem  with  the  scientific  tool/GUI 

https://www.dropbox.com/s/rar5h441umjdzgx/Appendixl.avi?dl=0 

Appendix  2:  Movie  demonstrating  solving  the  inverse  problem  with  the  scientific  tool/GUI 

https://www.dropbox.com/s/9qxvbul8naj6tpr/Appendix2.avi?dl=0 

Appendix  3:  Source  code  and  cross-platform  executable  of  the  scientific  tool/GUI 

https://www.dropbox.com/s/qwy5clzooyxthu5/Appendix3.zip?dl=0 
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